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We study the effects of spin-antisymmetric interactions on the stability of a Landau- Fermi liquid on 
the square lattice, using the generalized Pomeranchuk method for two-dimensional lattice systems. 
In particular, we analyze interactions that could induce instabilities of the so called spin-split type, 
that is when spin-up and spin-down Fermi surfaces are displaced with respect to each other. The 
phase space is studied as a function of the strength of the interaction V, the electron chemical 
potential \i and an external magnetic field h. We find that such interactions produce in general 
an enhancement of the instability region of the Landau- Fermi liquid. More interestingly, in certain 
regions of the V-fi phase space, we find a reentrant behaviour as a function of the magnetic field /i, 
similar to that found in recent experiments, e.g. in URu2Si2 and Sr3Ru207- 
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I. INTRODUCTION 

The low energy physics of weakly interacting fermions 
in three dimensions is mostly described by the Landau 
theory of the Fermi Liquid which is based on the exis- 
tence of single particle fermionic excitations, or quasipar- 
ticles, with a long lifetime at very low energies^— . How- 
ever, moving away from three dimensions, the situation 
changes drastically: in one dimensional systems Landau's 
quasiparticles are typically unstable, giving rise to the so 
called Luttinger liquid^. On the other hand, two dimen- 
sional lattice models are far more complicated to treat, 
since conventional perturbation theory is not applicable, 
and many instabilities can occur 5-26 . 

The appearance of broken symmetry phases due 
to electron interactions within the Landau description 
of a Fermi liquid has been the subject of numerous 
investigation s^ 16 ! 20 ! 27 ! 28 . In particular, Pomeranchuk in- 
stabilities leading to nematic phases have been analized 
by Yamase and Kohno in Ref. [25|. Of special interest 
are those phases that are still described by Landau's the- 
ory, such as metallic ferromagnetism or incommensurate 
density waves. In ref. 29 a new type of instability was 
discussed that lead to a phase in which the two Fermi 
surfaces for up and down spins get a relative shift in 
momentum space. This so-called spin-split phase (or a- 
phase 20 ), which breaks reflection symmetry but preserves 
time reversal (as opposed to a ferromagnetic metallic in- 
stability) was shown to appear in the presence of inter- 
actions that arise naturally in a tight binding description 
of interacting electrons. The same author have analyzed 
the competition of this kind of interactions with Rashba 
and Dresselhaus ones 30 . 

The subject has regained interest more recently, due to 
its possible relevance in the description of different phase 
transitions observed experimentally, such as the so-called 
hidden order transition in URu2Si2. More precisely, in 



ref. [3l| the possible relevance of spin-antisymmetric in- 
teractions in the stability of the Fermi liquid and the 
relation with this hidden order transition has been dis- 
cussed. Recent measurements, showing rotational sym- 
metry breaking in the hidden-order phase 3 - 2 -, motivated 
the proposal of a nematic emergent phase 3 - 3 -, in contrast 
to previous scenarios 3 -^. 

It is the purpose of this paper to analyze further this 
issue by applying a method that has been recently de- 
veloped in refs. Hi], HH, S, which is specially suited to 
study the Fermi liquid instabilities in two-dimensional 
lattice systems. This technique allows to analyze the 
interplay between an external magnetic field and the un- 
derlying lattice. In particular, we apply the extension of 
the above mentioned method to the zero-temperature 3 -^ 
and spinful case 3 - 6 - to the class of interactions studied in 

EH. 

Our results show two important characteristic features. 
On the one hand Pomeranchuk instabilities appear when 
the Fermi surface deformations are of the spin-split type 
and this leads in general to an enhancement of the insta- 
bility region. On the other hand, and more importantly, 
with this kind of deformations a reentrant behaviour in 
a magnetic field is observed, similar to what is observed 
in e.g. URu2Si2 and Sr3Ru207. These results should be 
compared with the analysis presented in |2l| where reen- 
trant behaviour was observed analyzing nematic defor- 
mations of the Fermi surface. Here we obtain a similar 
behavior, without imposing nematic deformations, for a 
different kind (spin-split) of Fermi surface deformations. 

The paper is organized as follows: in Section [IT] we 
present the Hamiltonian of the system and its descrip- 
tion as a Landau Fermi liquid, in Section [III] we briefly 
review the generalized Pomeranchuk method for lattice 
systems 3 -! — , and in Section [IV] we present the main re- 
sults of this paper. Finally in Section [V] we discuss the 
conclusions and outlook. 
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II. THE MODEL 

We study a system of fermions in a two-dimensional 
square lattice with a spin-antisymmetric interaction^. 
The Hamiltonian is given by 

H=^£(k)cl a Ck a + ^V^2 Jkk'(4a^^)(Ck'a' <7 « , i8 , ^ / ) ' 
k,o; k,k',i 

(2.1) 

where a % are the Pauli matrices and a = ±1/2 denotes 
the spin z-projection. We consider in this paper the dis- 
persion relation e(k) which corresponds to nearest neigh- 
bor hopping 

e(k) = -2t(cos k x + cos k y ) , (2.2) 

the analysis of more general cases being straightforward. 

For the interaction, following^! we choose J kk / in the 
antisymmetric spin channels 

(cos6> kk /). (2.3) 

i 

where # kk / is the angle between k and k'. For the sake 
of simplicity we consider in this paper a single /-channel, 
which we fix to I = 1, but our analysis can be easily 
generalized to higher /-channels. 

We take into account the presence of an external mag- 
netic field h by adding to (|2.ip the Zeeman interaction 

H mag = 2h^a c ka c ka • (2.4) 

k,Q! 

To carry the analysis further, we assume that the above 
system has a Fermi liquid phase that can be described 
using Landau's theory. Then we can write the change in 
the Landau free energy £7, associated with a change (5n kQ , 
in the quasi-particles distribution function n kcn as 




(2.5) 

where we have defined a spin dependent chemical poten- 
tial by reabsorbing the external magnetic field 

Ha = M + (2.6) 



Within the mean field approximation 5 , the dispersion re- 
lation in (|2.5p can be assumed to take the same form as in 
(|2.2p with renormalized coefficients. On the other hand, 
the "interaction function" f aa r(k, k') takes the form 

f aa , (k, k') = 2KJ kk , (2<W - 1) (2.7) 

The hypothesis of existence and stability of a Fermi liq- 
uid phase can now be tested: an instability is diagnosed 
whenever an excitation #n ka leads to a negative Landau 
free energy (|2.5|) . The way of parameterizing the space 
of excitations and exploring it in the search for possible 
instabilities is the core of the generalized Pomeranchuk 
method, and will be sketched in the next Section. 



III. THE METHOD 

In this section we give a brief description of the gen- 
eralization of Pomeranchuk's method to lattice systems, 
that was originally proposed in ref. 34 and further devel- 
oped in [HH and [36|. Here we concentrate in the spin 1/2 
case at zero temperature, with a dispersion relation (|2.2j) 
and interaction function (|2.7p and we refer the reader to 
ref. [36| for details. 

The ground state of a spin 1/2 system at zero temper- 
ature is described by its Fermi surfaces. With the help of 
the function # a (k) = fi a — s(k), the Fermi surface for each 
spin component can be defined as the loci in momentum 
space of the points at which # a (k) vanishes 

0a(k) = O. (3.1) 

The occupation number n kQ! vanishes outside the corre- 
sponding Fermi surface and takes the value 1 inside. 

At zero temperature, an excitation of the system £n ka 
takes the values 0, ±1 at different k points, and in con- 
sequence it can be characterized by a deformed Fermi 
surface. The latter can be represented by a perturbed 
version of equation (j3.ip 

g a (k)+5g a {k)=0 (3.2) 

In order to find a useful parameterization for 5<jf a (k), we 
change variables from (k x ,k y ) to (g a ,Sa), where s a = 
s a (k) is defined as orthogonal to # a (k), i.e. tangential 
to the corresponding Fermi surface, and performing a 
complete turn around it while running on the interval 
[— 7r,7r). The associated Jacobian reads 

Now if we choose 

Sg a (k) =^l^|^ 1 ^+^ 2 ^^J /3 j , (3.4) 

with e a p the Levy-Civita antisymmetric symbol en = 
€22 = 0, ei2 = —621 = 1, implying that the factor that 
multiplies i(j 2 is perpendicular to J a in a space. We have 
the deformations of the unperturbed Fermi surface corre- 
sponding to the a component (g a = 0) parameterized in 
terms of two arbitrary functions V ;1 ( 5 «) an d ^ 2 (s a ) of a 
single variable s a . The generality of this decomposition 
can be understood by replacing it into the expressions for 
the variations of the magnetization and the total number 
of particles. When ^(s) has no constant term it parame- 
terizes a simultaneous deformation of both Fermi surfaces 
that do not change the magnetization of the system nor 
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the total number of particles. Stoner instabilities, due 
to fluctuations in the number of particles, can be also 
included by allowing a constant term in ^(s). On the 
other hand ip 2 (s) takes into account deformations that 
include spin flips, changing the total magnetization. 



Inserting the above decomposition into the variation of 
the Landau free energy (j2.5p . after some straightforward 



algebra**, we get 



SQ = Jds Jds' ^ a (s)K ab (s,s f )^ b (s f ) 



(3.5) 



a, 6=1, 2 



As we have chosen both variables s a as running in the 
same interval, we are allowed to omit the index a in the 
dumb integration variables. Here K ab (s,s') is a 2 x 2 
matrix kernel given by 



K 11 (s, s') ■ 



K 22 {s,s f ) ■ 



K 12 {s,s')- 



J2 € ^ij aJp J n S{a - a') + ^e^J a (s)J p {s)J- 1 (s')J s (s')f a ds,s') 



2EaJ2WE a J!('') 



a, /3,7 



a, /3,7, 5 



- (j2e^J 2 J p 5(s-s') + E e^J a (s)J e (s)J y (s') 2 fa-r(s,s') 



2E a J|(«)Ea^(«') 



ot,/3 



ot,(3 ,-y 



(3.6) 



To find the instabilities of the Fermi liquid, we consider 
its energy ([3. 5ft as a quadratic form acting on 

the 2-vector function ^ a (s) that parameterize the de- 
formations, and then search for a choice of ip a (s) that 
turns it into negative values. A simple way to do that 
is to apply the Gram- Schmidt orthogonalization proce- 
dure to an arbitrary starting basis {^n( 5 )}nGiN, to obtain 
a new basis of orthogonal functions {£n( s )}neM satisfy- 
ing (£ n |£m) (X- ^nm- If at a given step of the procedure 
a 2-vector function £f (s) is found whose pseudonorm is 
negative, we say that we have found an instability in the 
i-th channel. In our case, we choose the initial basis as 



cos( 



ns > 
2 ) 







sin(^) 






cos(f) 



(3.7) 



if n even and a = 1 



if n odd and a = 1 



if n even and a = 2 



if n odd and a = 2 



IV. RESULTS 

By exploring the set of possible Fermi surface defor- 
mations as described in the previous section we obtained 
the unstable regions in the phase space. 




I 

In Fig. [1] we show deformed Fermi surfaces for both 
spin projections, corresponding to the first modes that 
lead to instabilities. As our method works with in- 
finitesimal deformations, the surfaces shown in the figure 
are schematic, since the deformation amplitude has been 
augmented for the sake of illustration. 





FIG. 1. (Color online) Fermi surface deformations leading 
to instabilities of the Fermi liquid. In blue (red) we plotted 
the spin-up (spin-down) perturbed Fermi surfaces, the dashed 
lines corresponding to the unperturbed ones. The first modes 
for n/t = —2, V/t = —2 are shown, the upper row corre- 
sponding to h/t = 0, the lower one to h/t = 0.5. Note the 
relative displacement of the spin up and spin down Fermi sur- 
faces in the first row, which indicates that the instability will 
probably lead to a spin-split state. 
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FIG. 2. (Color online) Unstable regions (shaded) for external 
magnetic field going from h/t = to h/t = 3, in figs, (a) to 
(d). 
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FIG. 3. (Color online) Instability region as a function of ap- 
plied magnetic field; by increasing it, is possible to leave (red 
line) and re-enter (blue line) to the stable Fermi liquid phase. 



We explored an ample range of parameters /i/t, V/t 
and h/t and obtained the instability regions shown in 



Fig. [2](a)to(d). One can observe that the instability be- 
comes increasingly significant for fi = where one of 
the Fermi surfaces gets close to the van Hove singularity. 

A remarkable feature of the present phase diagrams is 
that, for some values of \i and interaction strength V, the 
Fermi liquid description breaks down only for intermedi- 
ate values of h. For instance, as shown in Fig. [3l if the 
system is filled up to fi/t = 2 and the interaction strength 
is V/t = — 0.23, then we have a stable Fermi liquid up to 
h/t ~ 1.4 and from h/t « 2.6 but such Fermi liquid de- 
scription does not hold when h/t takes values in between 
those limits. This means that, under such conditions, it 
is possible for the system to leave and re-enter the Fermi 
liquid phase, as the external field is increased. 



V. CONCLUSIONS 

In this paper we applied the method developed invi - 
to a lattice model with spin degrees of freedom and in the 
presence of interactions that break the spin symmetry^. 
The phase space was studied as a function of the interac- 
tion strength V, the electron chemical potential \i and an 
external magnetic field h. Our results show that Pomer- 
anchuk instabilities appear when the Fermi surface defor- 
mations are of the spin-split type (that is when spin-up 
and spin-down Fermi surfaces are displaced with respect 
each other in momentum space), confirming the studies 
in [29| and [3l|. This leads in general to an enhancement 
of the instability region. On the other hand, and more 
importantly, we observe a reentrant behaviour in a mag- 
netic field, similar to what is observed in e.g. URu2Si2 
and Sr 3 Ru 2 07 (see e.g. |22[ ). 
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